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Abstract. The brightest events in a time series of cosmological transients obey an observation 
time dependence which is often overlooked. This dependence can be exploited to probe the global 
properties of electromagnetic and gravitational wave transients [1, 2]. We describe a new relation 
based on a peak flux-observation time distribution and show that it is invariant to the luminosity 
distribution of the sources [3] . Applying this relation, in combination with a new data analysis filter, 
to Swift gamma-ray burst data, we demonstrate that it can constrain their rate density. 
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INTRODUCTION 



Transient astronomy is a rapidly emerging field which will benefit enormously as coor- 
dinated observations of sources emitting electromagnetic photons, gravitational waves, 
cosmic rays or neutrinos become commonplace. Transient sources, with peak emission 
durations much less than typical observational periods, will include events such as su- 
pernovae, gamma-ray bursts (GRBs) and gravitational wave emissions from compact 
Lj_ ! object inspirals [1]. 

^ For cosmological populations of sources, the brightness distribution is conventionally 

q ■ used to constrain their luminosity function (LF), their evolution in density [4] and for 

transient sources, their rate density [5]. Estimates of global properties are obtained by 
fitting the number - brightness distribution to models that include luminosity, source 
density and evolution effects. For the case of transient events an additional parameter is 
available - the event arrival times. 

The brightest events of an astrophysical transient population obey a simple, not well 
known, observation-time dependence. Figure 1(A) illustrates this dependence by way 
of the observation time - peak flux distribution of long duration (> 2 s GRBs (LGRBs) 
as observed by the Swift satellite 1 . Although GRBs are not standard candles and ob- 
servations of these transients suffer from selection effects, two things are immediately 
apparent from this plot: Firstly, the initial observations are of the intermediate energy 
sources - these events occur frequently and correspond to the high-z peak of the event 
distribution. Secondly, as observation time increases, rarer events from both z extremes 
of the distribution come into play, producing a triangular distribution in log-log space. 



1 http://swift.gsfc.nasa.gov/docs/swift/swiftsc.html 



The events that are missed through observational bias or background detector noise will 
not generally be amongst the brightest events in a distribution. Therefore, to extract in- 
formation about the underlying source population it could be beneficial to use the the 
brightest events in a survey. 

The temporal distribution of transient astrophysical populations of events has been 
described by the 'probability event horizon' (PEH) concept of [2], This method estab- 
lishes a temporal dependence by noting the occurrences of successively brighter events 
in a time series - shown by the highlighted events in Figure 1(A). By utilizing the fact 
that the rarest events will preferentially occur after the longest observational periods, 
it produces a data set with a unique statistical signature. To extract this signature from 
a body of data, Howell et al. [1], developed a PEH filter, which used in unison with a 
brightness - observation time relation, could be applied to obtain estimates of the rate 
density of events, tq, from a standard candle population of sources. In its simplest form, 
the PEH filter simply extracts successively brighter events as a function of observation 
time. However, for the case of GRBs and many other transient sources, there is a distri- 
bution in luminosities. This imposes a severe restriction on the technique. To overcome 
this problem, Howell et al. [3] derived a new peak flux - observation time relation that 
can be applied to distributions of transient astrophysical events and is independent of the 
form of the LF. This relation will be described in the next section. 



THE PEAK FLUX - OBSERVATION TIME RELATION 

Consider a distribution of events distributed in a Euclidean space by an event rate density 
tq and a LF <j>(L) (L m in < L < L max ). The observed peak flux, P (photons cm s ), or 
'brightness', distribution of events over an observation time T is a convolution of the 
radial distribution of the sources and their LF. For tq and (L) independent of position, 
setting the upper limit of distance r is the maximum range for which an event with 
luminosity L, produces a peak flux P. Integrating over the radial distance yields: 

N(>P)= T J°^ P -ift f~ hdlWjl, (i) 

3 V47T JL mm 

where the average solid angle covered on the sky has been accounted for by A£2/4tt. 
This is the familiar log AMog P relation, N(> P) <* P~ 3 > 2 , a power law independent of 
the form of the LF [6] . 

To see how the events are distributed in time, we note that, as the individual events 
will follow a Poisson distribution in time the temporal separation between events follow 
an exponential distribution, defined by a mean event rate R(r) = ro(A/3)nr 3 for events 
out to r. The probability for at least one event > P to occur in a volume bounded by r 
during an observation time T at constant probability £ is given by: 

&>(n>\\R{r),T) = \-tP { ri T = e. (2) 

For this equation to remain satisfied with increasing observation time: 

N(>P)=R(r)T = \\n(l-e)\. (3) 

Equations (1) and (3) for N(> P) combine to give the relation for the evolution of 
brightness as a function of observation time [3]: 
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FIGURE 1. (A) The observation time - peak flux distribution of long gamma-ray bursts as observed by 
the Swift satellite. Successively brighter events that have occurred during this period are highlighted. 
(B) PEH fits to the Swift LGRB data sample. The estimates are consistent with estimates by [7] of 
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This relation shows that for a Euclidean event distribution, a log P-log T distribution 
will have a slope of 2/3, independent of the form of the LF. We can use the log P-log T 
relation to produce curves defining the probability, £, of obtaining some value of peak 
flux, P, within an observation-time, T, for a given tq and (j>(L). 



A TIME-BASED SIGNAL PROCESSING FILTER 

Figure 1(A) shows that in its simplest form, the PEH filtered output consists of only a 
small sub-sample of the overall population. Thus, estimates of parameters obtained by 
fitting to a model PEH curve will inevitably suffer from low resolution. To increase the 
data sample, one can split a time series into smaller time segments, apply the filter to 
each and recombine the output [1]. The use of this technique is subtle, as splitting can 
increase the number of outlying events and can diminish the statistical signature of the 
output. In addition, there is a significant probability of a bright event occurring early in 
an observational period, further reducing the sample size. 

To increase the amount of usable data and to enhance the signature of the filter 
output, two physical principles can be exploited [3]. Firstly, the temporal cosmological 
principle, implies that for time scales that are short compared to the age of the Universe, 
there is nothing special about the time we make an observation. Thus, the PEH signature 
of a transient population of events is independent of when a detector is switched on. Time 
reversal invariance means that a time series of events can be treated as a closed loop. 
Therefore, the PEH filter can be applied to a data set from any start time and in both 
temporal directions. As the start time for the PEH analysis is now arbitrary, the filter can 
be applied in such a way that the brightest event can be set as the final event in a series, 
and the observational period will be the total length of the loop. These techniques can 
increase the PEH sample and significantly improve the statistical signature of the filter 
output. 



APPLICATION TO SWIFT DATA 

To demonstrate the log P- log T technique, Howell et al. [3] used the LGRB peak flux 
data recorded by the Swift satellite between 2004 December and 2007 April. Validation 
was achieved by fitting a 90% confidence log P- log T threshold to PEH filtered data, 
using recent estimates of source parameters from [7]. 

In this section, we explore the efficiency of the method as a practical tool for probing 
global source parameters, using a least-squares fit to the same data. Estimates of the 
parameters are compared with those of [7], who fitted to the brightness distribution to 
estimate the rate density and LF parameters of the LGRB population. 

To allow for a cosmological distribution of sources, equation (4) must be modified 
to allow for cosmic evolution. A standard Friedman cosmology is used to define a 
differential event rate, dR(z), in the redshift shell z to z + dz. The luminosity and flux 
will be related through z by a luminosity distance d^(z) [see, e.g., 2, 8]. In this case, 
solving equation (3) numerically, with P = L/47td^(z), will yield the cosmological 
log P- log T relation. We apply this fit to the same GRB data, setting £ to 50%, the 
mean value obtained from 10,000 simulated data sets. 

Figure 1(B) shows the log P- log T fit to the Swift LGRB sample, yielding tq = 
0.10± 0.02 Gpc~ 3 yr _1 at a l-<7 confidence level - this is within the range of tq = 
O.lO^QQg Gpc~ 3 yr~ ! obtained by [7] using four free parameters. The greater resolution 
of the log P- log T fit is a result of using only one free parameter, r$. 

The log P- log T fit utilizes the brightest events in a data stream, which on average 
correspond to the closest events [9]. Therefore, an accurate representation of the high-z 
selection function may not be essential. Additionally, the technique may be a useful 
probe of the high end of the LF. Preliminary studies [10] suggest the method may 
shed some light on the narrow observed distribution of BATSE peak energies [11, 12]. 
Finally, we note that the time dependence conveniently allows the method to be used as 
a predictive tool, to determine the probability of an event with some predefined energy 
to occur within some observation time [3, 13]. 
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